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When calls offered to a primary group of trunks find all of them 
busy, provisions are often made for these calls to overflow to other 
groups of trunks. Such traffic overflow systems have been of interest 
for a long time, but recently overflow systems that allow for some 
calls to be queued have been of importance. In this paper we analyze 
a traffic overflow system with queuing. The system consists of two 
groups, a primary and a secondary. We consider two cases which 
differ in the treatment of demands waiting in the primary queue. We 
present an analytical approach which is suitable if the secondary 
queue is large, or even infinite, and we contrast it with an earlier 
approach of ours which is more suitable if the secondary queue is not 
large. The analysis considerably reduces the dimensions of the prob- 
lem and simplifies the calculation of various steady-state quantities 
of interest. Our results include expressions for the loss probabilities, 
the average waiting times in the queues, and the average number of 
demands in service in each group. 



I. INTRODUCTION 

In an earlier paper, 1 we analyzed a traffic overflow system with 
queuing. The system consists of two groups, a primary and a secondary, 
with n k servers and q k waiting spaces, which receive demands from 
independent Poisson sources S k with arrival rates X* > 0, k = 1 and 2, 
respectively, as depicted in Fig. 1. The service times of the demands 
are independent and exponentially distributed with mean service rate 
H > 0. If all no servers in the secondary are busy when a demand from 
S2 arrives, the demand is queued if one of the q 2 waiting spaces is 
available; otherwise, it is lost (blocked and cleared from the system). 
Demands waiting in the secondary queue enter service (in some 
prescribed order) as servers in the secondary become free. 
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Fig. 1— Mean flow rates for an overflow system with queuing. 

If all /ii servers in the primary are busy when a demand from Si 
arrives and there is a free server in the secondary and no demands 
waiting in the secondary queue, the demand is served in the secondary. 
If there are no free servers, then the demand is queued in the primary 
queue, if one of the qi waiting spaces is available; otherwise, it is lost. 
Two different cases were considered for the treatment of demands that 
are waiting in the primary queue. In case I, a demand waiting in the 
primary queue may enter service either in the primary, when a server 
becomes free, or in the secondary, if a server becomes free and no 
demands are waiting in the secondary queue. In case II, no overflow is 
permitted from the primary queue, so that a demand in the primary 
queue must wait for a server in the primary to become free. 

In this paper, we use an alternate approach to analyze the same 
overflow system. This alternate approach is useful, in particular, when 
qi is large, or even infinite. Which approach is preferable depends, 
essentially, on the relative size of qz to rii . Other approaches to the 
problem were discussed in the earlier paper. 1 Letp,> denote the steady- 
state probability that there are i demands in the primary and j 
demands in the secondary, either in service or waiting. These proba- 
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bilities satisfy a set of generalized birth-and-death equations, which 
take the form of partial difference equations connecting nearest neigh- 
boring states. Here we carry out an analysis that reduces the dimen- 
sions of the problem, which may be considerable in cases of interest. 
An analogous reduction was obtained in the earlier analysis. 1 

The basic technique is to separate variables in regions away from 
certain boundaries of the state space, the elements of which are {i, j). 
These regions are depicted in Fig. 2a. The analogous regions corre- 
sponding to the analysis of the earlier paper 1 are depicted in Fig. 2b, 
for comparison. The separation of variables leads to two sets of 
eigenvalue problems for the separation constant. The eigenvalues are 
roots of polynomial equations. The probabilities p,> are then repre- 
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Fig. 2— Boundaries of regions in state space for the analysis of (a) this paper, and (b) 
the earlier paper. 
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sented in terms of the corresponding eigenfunctions. The constant 
coefficients in these representations are determined from the boundary 
conditions and the normalization condition that the sum of the prob- 
abilities is unity. In general, these constants have to be determined 
numerically. In case II, additional accessible states are in the region 
bordered by the broken lines in Fig. 2, since in this case demands wait 
in the primary queue even if a server is free in the secondary. In this 
additional region, the probabilities />,, are expressed, as in the earlier 
analysis, 1 in terms of a fundamental solution of a partial difference 
equation. 

There are various steady-state quantities of interest, which may be 
expressed in terms of the probabilities />,> . The quantities include the 
loss (or blocking) probabilities, the average waiting times in the queues, 
and the average number of demands in service in each group. These 
quantities may be expressed directly in terms of the constant coeffi- 
cients which occur in the representations for the probabilities p,y. Thus 
the steady-state quantities of interest may be calculated directly, once 
the coefficients have been determined from the boundary and normal- 
ization conditions, without having to calculate the probabilities p,y. 
Here again the reduction in the dimensions of the problem is valuable. 

We first consider case I. The representation of the probabilities /?,> 
in terms of the eigenfunctions is discussed in Section II, and the 
boundary and normalization conditions are considered in Section III. 
Various steady-state quantities of interest are calculated in Section IV. 
Next, we consider case II and discuss the representation of the prob- 
abilities pij and the boundary and normalization conditions in Section 
V. The corresponding steady-state quantities of interest are considered 
in Section VI. Properties of the eigenfunctions which occur in the 
representations of the probabilities p,y are given in the appendix. 

II. REPRESENTATION OF SOLUTION: CASE I 

Let p^ denote the steady-state probability that there are i demands 
in the primary and j demands in the secondary, efther in service or 
waiting. These probabilities satisfy a set of generalized birth-and-death 
equations, 2 which may be derived in a straightforward manner. We 
define the traffic intensities 

ai = Ai/ju, a 2 = A 2 //i, (1) 

and let 

ki=n! + q lf k 2 = n 2 + q2. (2) 

It is convenient to introduce the function 



fl, t > 
10, € < 



x-h;/ ; ^ J (3) 
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as well as the Kronecker delta 

Ht'^J (4> 

Then, 1 for case I, it is found that 
[ai(l - Si^xj-ru,) + 02(1 - fi/Aj + min(i, /ii) + min(y, n 2 )]pv 
= (1 - Xhi-*iXH-i-j)l a *0- ~ s io)Pi-ij 

+ (1 - 5,*,)min0' + 1, n 2 )pu+i] 

+ (1 - 5>o)[aifiw 1 x« 2 -> + 02(1 - X'-i-n,Xi2->)]P'.y-i. 

+ (1 - 6iA,)[(l - x»-»,>&4-i-y)m"i(i + 1. »i) + ^x.-nAnJP'+w. (5) 

for < 1 < k\, </' < k 2 . These equations were constructed to imply 
that 

Pa = 0, m + 1 < i < k u <y < n 2 - 1, (6) 

since it is impossible for demands to be waiting in the primary queue 
when a server is free in the secondary. The normalization condition is 

1 £ Pii+1 ZW-I. (7) 

i-O j=0 i-0 j=ri2 

We assume that q 2 > 1. Then, for < i < &i and n 2 + 1 <y < £ 2 , the 
variables in (5) may be separated, and there are solutions of the form 
otifij where 

[ai(l - 8ik t ) + min(i, nt) + fc]a, 

= ai(l - 8io)oti-i + (1 - 6iA,)min(i + 1, ni)a,-+i (8) 

for < i < /si, and 

[a 2 (l - 8 Jk2 ) + n 2 - te]0j = a 2 pj-i + n 2 (l - S/kJfai (9) 

for n 2 + 1 <y < A> 2 , and k is a separation constant. The solution of (9) 
may be expressed in terms of Chebyshev polynomials of the second 
kind, 3 U/(x). It is convenient to define 

w= «';fe), (10) 

\«2/ V 2Va 2 n 2 / 
and 

fcOt) = Vkr-jiK) ~ *k 2 -J-l(K). (11) 

The properties of these functions which we will need are given in the 
appendix. We note here, however, that Uo(x) = 1, U-\(x) = and 
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4>k 2 M = 1. It follows from (9) and (102) that /?> is proportional to </>,(«) 
for 7i 2 < y < k 2 . 
For < i < Hi - 1, (8) implies that 

(ai + i + K)ai = ai(l - 5io)a«-i + (i + Da l+ i. (12) 

The solution of (12) may be expressed in terms of Poisson-Charlier 
polynomials. 4,5 We here denote the solution of (12) for which a = 1 by 
s,(k, ai). The properties of s,(k, a) which we will need are given in the 
appendix. For /ii < i < ki , (8) implies that 

[a,(l - 8a t ) +ni+ /c]a, = dica-i + #ii(l - 5,A,)a,+i. (13) 

Corresponding to (10) and (11), we define 

QAk) = (*Yu J ai+ £l*) (14) 

and 

0,0c) = GftH<*) ~ Qft,-f-iW. (15) 

It follows from (13) and (111) that a, is proportional to 0,(k) for /ii — 1 
< t < k\ , and we note that 0*, (/c) = 1. 
Consequently, we take 

f *•■(«, ai)^„,(K), 0<i</ii, . . 

a ' \ «„,(«, ai)$i(K), m-l<t<A>,, u °' 

where 

s n ,_i (k, ai )0 n , (k) = s„j (k, ai )0„,-i (k). (17) 

With the help of (15), (96), (97), and (110), this equation may be 
written in the form 

k[s„,(1 + k, ai)fl 9l (/c) - s„,-i(l + k, ai)fl,,_i(ic)] = 0. (18) 

The expression in the square brackets in (18) is a polynomial in k of 
degree k\ = n.\ + q\. It was shown 1 that its zeros are negative and 
distinct, and we denote them by K r , r = 1, • • •, ki. We also let *o = 0. 
It was also shown 1 that the zeros of ^ M are positive. Hence, 
$n,(jCr) f* 0, r « 0, • • •, k\. Moreover, it follows from (118) that 



lim +' M 



9z-»°° </>/i 2 (fr) 

for fixed j, where 



"2( €r+v /^TT)_ 



r-1, ..-,*,, (19) 



er= (a 2 + n i -0 > (2Q) 

2va2/i2 
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and the positive square root of el - 1 is taken in (19). We note, from 
(109), that 

I 10,9 

e r + Ve? - 1 > -v /— , r=l, ...,*i, (21) 

since K r < 0, r = 1, • • • , k\ . On the other hand, from (105), 

*d®-=hV>* (22) 

and, as expected, the condition a 2 < /i 2 is necessary for a steady-state 
solution when q 2 = °o. 

In view of the above results, we may represent the probabilities p, y , 
for n 2 &j "& k 2 , in the form 

£ b r Si(K r , ai)s„ 2 (-K r , a 2 )0 nj (K r ) / / , < t < Hi, 

r-0 <P/i 2 (' C '-) 

^ _ *' <b (k ) 

£ 6 r s„,(/c r , ai)s^(-K r , a 2 )di(K r ) / /. , RiStsAi, (23) 

where the constants 6 r are to be determined. The reason for the factor 
s^—Kr, a 2 ) will be apparent shortly. We note, from (94), that 
Sn 2 (—K r , a 2 ) > 0, since K r < 0, r = 0, •■•,&. 

For < i < Tii — 1 and < j < rc 2 — 1, the variables in (5) may be 
separated, and there are solutions of the form y&j where 

(ai + i + ri)yi^ ai(l -8i Q )y i - 1 + (i + l)y i+u OSiSiii-1, (24) 

and 

(oa +> - Tj)5 y = a 2 (l - fi^)«,-i + U + W«, <y < n 2 - 1, (25) 

and 7j is a separation constant. It follows from (93) that y, is propor- 
tional to s,(tj, ai) for < i < rii and 5, is proportional to Sj(—q, a 2 ) for 
< j < n 2 . Now s„ 2 (— 7j, a 2 ) is a polynomial in tj of degree n 2 . It was 
shown 1 that its zeros are positive and distinct, and we denote them by 
i\(, «f = 1, • • •, n 2 . Then 

fe,(-i|A as) - 0, <?= 1, • • -, /i 2 . (26) 

Consequently, we represent p, y , for < i < ni and < / < n 2 , in the 
form 

/?«>= £ 6 r s,(fc r , ai)s>(-ic r> a 2 )0 ni Oc r ) + £ c^s,(rj A ai)s,(-7j,, a 2 ), (27) 

where the constants C/ are also to be determined. Note that the 
representations in (23) and (27) agree fory = n 2 , < i < m, in view of 
(26). 
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It remains to satisfy the boundary conditions at i = ri\ , < y < 
/12 — 1 and aty = n 2 , < i < k\, as well as the normalization condition 
(7). This is done in the next section. 

III. BOUNDARY CONDITIONS: CASE I 

From (5), the boundary conditions at i = n\ imply that 

(ai + a 2 + ni +j)Pn u 

= aipn x -\j + (ai + a 2 )(l - Sjo)p nj j-i + (j + l)Pn lv/+ i, (28) 

for <y < n 2 — 1. Also, the boundary conditions aty = n 2 imply that 

(d + a 2 + i + n 2 )pw 2 , 

= ai(l - 8io)p i -i,n 2 + tHPi^-t + (i + l)p,+i^2 + n 2 p,- l n 2+ i, (29) 

for < i < ni - 1, 

[ai(l — 5 9l , ) + a 2 + ni + «2]pn,,n 2 = aiPm-i.il, 

+ (ai + a 2 )Pn,,n 2 -l + (/ll + ^2)(1 - 5 9 „o)Pn 1 +l,n a + n 2 Pn 1 jt^+l, (30) 

and, if ^i > 1, 

[ai(l - &*,) + a 2 + ni + n 2 ]p«> 2 

= aipi-i^ + (/ii + n 2 )(l — &*, )/>i+i,n2 + nzPbig+ii (31) 

for /ii + 1 < i < Ai . 

If we substitute (27) into (28), we find, after reduction with the help 
of the recurrence relations in the appendix, that 

£ brdn^Kr^KrSn^l + Kr, Oi )«/(-«,-, 2 ) + aifi„, (fCr^l )S,(~1 - Kr, a 2 )] 
r-0 

"2 

+ £ CtffaAi,(l + "m* a\)sj{—i)t, a 2 ) 

+ ai**, (i?ft «i)*/(-l - fe a 2 )] = 0, (32) 

for <y < n 2 - 1. We remark that the first sum in (32) may be written 
in a different form with the help of the relationship 

K r [Sn,(l + K r , ai)0„,(lCr) + Sn^r, ai)J2 gi -i(jC r )] = 0, (33) 

which follows from (15), (97), and the fact that K r , r = 0, • • •, ki, are 
the roots of (18). 
If we make use of (23) and (27) in (29), we find after reduction that 



a 2 X b r Si(K r , a\)8 n ,(Kr) 



Sn. z (-K r , a 2 ) Sn^-li-Kr, d 2 ) 

VnzKKr) 

«2 

- a 2 X c f Si(f]/!, a\)Sn 2 -i(— t\tt a 2 ) = (34) 
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for < i < n\ — 1. The first sum in (34) may be written in a different 
form with the help of the relationship 

a 2 [Sn. 2 (-K r , a 2 )<t>n 2 -\(Kr) ~ 8^-\(—K ri Cfe^fcr)] 

= -K r [Sn 2 (l ~ K r , a 2 )% 2 (K r ) ~ ft*-l(l - K r , a 2 )% 2 -l far)], (35) 

which follows from (11) and the recurrence relations in the appendix. 
We note that k = 0. Also, for q\ > 1, if we make use of (23) in (31), it 
is found after reduction that 



r=0 



a 2 8i(K r ) 



<f>n 2 -l(K r ) 
<t>n 2 (K r ) 

- n 2 (l ~ 8 ikl )6 M (Kr) 



=0, (36) 



for ni + 1 < i < k\. 

If we substitute (23) and (27) into (30), and make use of (17), (102), 
and (111), we find that 



£ brSn^Kr, Cli)Sn. 2 {—Kr, Cl 2 ) 
r=0 



aiOn^Kr)- 



<t>no-\{*r) 



- tl 2 {\ ~ S qi fi)0 ni +l(K r ) 



fynMr) 

- (ai + a 2 ) X b r s„ x (Kr, aJSnt-A—Kr, a 2 )0 nx {K r ) 

r-0 

"2 

- (ai + a 2 ) £ CtSn^Tis, aiUnt-ii—n*, a 2 ) = 0. (37) 



This condition may be shown to be redundant, by summing (32) from 
/ = to 7i2 — 1, (34) from i = to n\ — 1 and, for q\ > 1, (36) from i = 
ti\ + 1 to ki, and adding. 

The constants b r , r = 0, • • • k\ and c (y *f= 1, • • •, n 2 are determined 
by (32), (34), and (36) only to within a multiplicative constant, which 
is determined from the normalization condition (7). But, from (23), 
with the help of (15), (18), (97), and (98), it is found that 

t Pv = b oSn2 (0, a 2 )^-[s ni a, ai)n,,(0) - s^-id, o,)O gi .,(0)], (38) 

for n 2 <j<k 2 . Then, from (7) and (27), with the help of (11) and (98), 
it follows that 

bos n .,(0, (feP^^d, oi)Q„(0) - s„,-i(l, 01)0^-1(0)] 

<P* 2 (U/ 

+ X 6 r S n ,(l + K r , ai)S^,-i(l - K r , a 2 )0n x {^r) 



+ J c,s„,(l + t,,, ai)s„ 2 -i(l - t,,, a 2 ) = 1. (39) 
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We note, from (104) and (105), that 

,. ¥,,(0) n 2 c ^ ,. m 

hm ' =- r, for a 2 < n 2 . (40) 

92 ^=o ^ 2 (0) (zi2 - a 2 ) 

Once the constants b r and c, have been determined, the steady-state 
probabilities p<, may be calculated from (23) and (27). We remark that 
the number of constants to be determined is only k\ + 1 + n 2 , whereas 
the number of probabilities p y is (q 2 + l)(fci + 1) + n 2 (ni + 1), which 
in general is considerably larger. It is of interest to compare the 
approach taken in this paper with that adopted in the earlier paper 1 
(see Fig. 2). There the variables in (5) were separated in the region 
< i < rii - 1 and < j < k 2 , leading to k 2 + 1 eigenvalues and 
eigenfunctions, and, for <?i > 1, in the region n { + I < i < ki and 
n 2 + 1 <y : < k 2 , leading to q\ eigenvalues with corresponding eigen- 
functions vanishing at i = n x . In the corresponding representation of 
the probabilities py in terms of the eigenfunctions, there are thus k 2 
+ 1 + <7i constants to be determined. In view of (2), which approach is 
preferable depends on the relative size of q 2 to n\. The approach 
adopted in this paper is suitable if q 2 is large, or even infinite. 



IV. SOME STEADY-STATE QUANTITIES: CASE I 

We proceed now to the calculation of various steady-state quantities 
of interest. These quantities are depicted in Fig. 1, which indicates the 
mean flow rates. The loss probabilities L\ and L 2 are given by 

Li - 2 P*u/. L 2 =t P i, k 2 , (41) 

and the probabilities that a demand from the primary, or secondary, 
source is queued on arrival are 

Q 1 m (1 - Vo) I i PU, <?2 - i j, Pij- (42) 

The probability that a demand arriving from the primary source 
overflows immediately is 

/12 = "S p ni j. (43) 

j-o 

Since the mean service rate is /t, the mean departure rate from the 
primary queue to the primary servers is 

*l *2 

Rn = »i/i(l - S 9l ,o) S S PV. (44) 

i=n| + l /— 14 

while the mean rate of overflow from the primary queue to the 
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secondary servers is 



k. 



Rl2 = #l 2 /l(l - S-7,.0) E Pi"2- ( 45 ) 



i=n,+l 



The mean departure rate from the secondary queue is 

*l *2 

R22 - /i 2 /i 2 £ p,>. (46) 

1=0 y=n. 2 +i 

It may be verified from (5) that 

Rn + R\2 = \lQl, i?22 = A 2 Q 2 . (47) 

These relationships hold since, in the steady state, the departure rates 
from the queues are equal to the arrival rates to them. 

The average number of demands in the primary and secondary 
queues are 

Vi= I I (i-ni)pij, V 2 =i t (j-n 2 ) PiJ . (48) 

'-"1 j-n-i i=0_/=/io + l 

Also, the average number of demands in service in the two groups are 

Xi = s I iPv + Mi - K.o) i i Pij < 49 > 

i-O y=0 i = n, + l>-/i2 

and 

x 2 = "i Yjpo + 12 i i Pii . (50) 

i-O y-o i'=0>-n. 2 

If we apply Little's theorem 2 to the primary and secondary queues, we 
find that the average waiting times of the demands which are queued 
in the primary or in the secondary are given by 

w^rrr^ 1 - 1 *' W2 = tt>-> (51) 

Ait^i A2^/2 

respectively, independently of the order of service within each queue. 
Also, if we apply Little's theorem to the primary and secondary groups 
of servers, we obtain 

A1U-L1-I12) -R i2 = nX 1 , \2(l-L 2 )+\J l 2 + Ri2 = nX2, (52) 

since the mean service rate is fi. 

The steady-state quantities of interest may be expressed in terms of 
the constants b r and c, with the help of the representations in (23) and 
(27). From (41), it is found, with the help of (11), since 0*,(/c) = 1, that 

L\ = X brS ni {K r , ai)S„.,(-K r , Og) "* T , (53) 

r-0 <Pn 2 KK r ) 
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and, with the help of (22) and (38), that 

L 2 = bo(^) s n2 (0, a 2 )[s ni (l, ai)Q 9l (0) - s ni -i(l, ai)Ofc-i(0)]. (54) 
From (42), it is found, with the help of (11) and (15), that 

Ql = S brSnMr, Oi^-*,, a*)[Q qi (K r ) ~ l]" 2 ^, (55) 

r-0 <Pn 2 {Kr) 

and, with the help of (38), that 

Q , [* g2 (Q) - 1] 
Q2_6 ° 0,(0) 

.3,(0, a 2 )[s n ,(l, ai)n 9l (0) - «b,-i(l, o,)fl, r i(0)]. (56) 

Next, from (43), with the help of (98), it is found that 

I\2 = S brS ni (K r , d^S^-lH - fC r , a 2 )6 ni (K r ) 
r=0 

"2 

+ 2 c,s ni (T)t, a,)s,_,(l - i\ t% a 2 ). (57) 
Also, from (44) and (45), if we make use of (11) and (15), we obtain 

*' ¥ (K r ) 

Rn = riifi £ brS ni (Kr, ai)s,(-K r , a 2 )tt Ql -i(>c r ) " 2 \ (58) 

and 

*' 
Ri 2 = n 2r i £ 6rS„,(K r , ai)s n2 (-fc r , a2)fi gi -i(/cr). (59) 

r=0 

From (46), with the help of (11), (38), (56), and (104), it follows that 
R 22 = X 2 Q 2 , as expected. The first relationship in (47) provides a useful 
numerical check, in view of the expressions in (55), (58), and (59). 
From (48), with the help of (11), (114), and (115), it is found that 



Vi = 6os,,,(0, ai)s,(0, a 2 )[ — A 9l — 



+ £ ~~ S nMr, ai)S,(— K r , d 2 ) 
r-1 K r 

• {a,[fi 9l ( Kr ) - 1] - niQ gi -iM) ^t^T, (60) 

<Pn 2 \Kr) 



where 



^"^ Wlff + D. €-t (61) 
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Also, with the help of (22) and (38), it is found that 

V 2 = 6oA 92 (^)s n2 (0, «)[*,(!, a,)O 7l (0) - s nr ,(l, aOO^-ifO)]. (62) 

It follows from (51), (56), and (62) that 

w _ A q2 (a2/n 2 )<t>n 2 (0) (6S . 

2 A 2 [* 92 (0) - 1] * 

Next, from (49), with the help of (7), (11), (98), and (99), it is found 
that 

n 2 

Xi = «i - X CA,-i(2 + r),, ai)s^,-i(l - i^, 02) 

*i 

- J 6 r S„,-l(2 + K r , ai)0 ni M 



% 2 (Kr) 
l^(/Cr) 



S n ,-l(l - Kr, a 2 ) + Sn,(-K r , d 2 ) 



. (64) 



Similarly, from (50), it is found that 



n 2 

X 2 = n 2 - £ <V*n,(l + T)/, ai)s^_i(2 - 7j,, a 2 ) 
e-\ 

*. 
- X &rSn,(l + Kr, ai)s^,-i(2 - K r , a 2 )0„, (/c r ). (65) 

r-0 

In view of (53), (54), (57), (59), (64), and (65), the relationships in (52) 
provide a useful numerical check. 
From (107) to (109), it follows that 

r *«* M (e r + ve?-l) . , fifi , 

hm — ^- — - = . , , r = 1, • • • , ki , (bb) 

9 2 ^co (|) n2 (jCr) ( €r + y/ € 2 r - 1 - Va 2 /n 2 ) 

where e r is given by (20). With the help of (40) and (66), we obtain the 
results for the limiting case q 2 = oo, with a 2 < n 2 , in (53), (55), (58), 
(60), and (64). From (54), as expected, L 2 — > as q 2 — > oo, for a 2 < n 2 . 
Also, since <J>* 2 (0) = 1, it follows from (22) and (40) that 

Um [ *'' (0 >7 1] - 7 -^- r . fora 2 < n2 , (67) 

,,-.« <J>„ 2 (0) (n 2 - a 2 ) 

and from (61) that 

lim A„ ? I — ) = — 2 — -2 , for a 2 <n 2 . (68) 

92 -.=o \n 2 ) (n 2 - a 2 ) 

From (67) and (68), we obtain the limiting results in (56) and (62), and 
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from (1) and (63), we have 

lim W 2 -—. -, for a 2 <n 2 . (69) 

«7 2 -.cc /Ll(/l2 — 0.2) 

We remark that the steady-state quantities of interest may be 
calculated directly, once the constants b r and c f have been determined, 
without having to calculate the probabilities p,y. 

V. REPRESENTATION AND BOUNDARY CONDITIONS: CASE II 

We now consider the second case, in which no overflow is permitted 
from the primary queue, so that a demand in the primary queue must 
wait for a server in the primary to become free. Since case II differs 
from case I only in the treatment of demands waiting in the primary 
queue, we assume that q\ > 1, as well as q 2 > 1. The steady-state 
probabilities p,> now satisfy the equations 1 

Oi(l - 5,a,X/-/. 2 ) + «2<1 - 8jk 2 ) + min(i, /ii) + min(y, n 2 )]p,> 

= cti(l - 6,o)(l - X'-«,-iX'»2-i->)/ > '-i.7 
+ (1 - 5 / o)(aiX«-/« 1 X'' 2 -7 + a2)pu-i 
+ (1 - 8ik,)mm(i + 1, ni)p i+i j 
+ (1 - 8jk 2 )min(j + 1, n 2 )p iJ+ i, (70) 

for < i < ki t —J — k 2 . The normalization condition is 

1 SP.y=l. (71) 

1=0 y=o 

Now (5) and (70) are identical for < i < /ii — 1, <y < n 2 and for 
< i < Ai, /12 + 1 <7 < ^2. It follows that the representations in (23) 
and (27) are still valid, although the constants b r and c, will, of course, 
be different. However, the boundary conditions (29), which correspond 
toy = n 2 for < i < /ii — 1, are still satisfied, and hence (34) still holds, 
for < i < n\ — 1. The remaining boundary conditions are different. 
Also, from (70), for rii + 1 < i < ki and <y < n 2 - 1, we have 

(ai + a 2 + rii +j)pij 

= (ai + a 2 )(l - 8j )pij-i + /ii(l - Sikjpt+ij + (j + l)pij +l . (72) 

We define 1 the quantities II my , for m,j = Q, 1, • • • , as the solutions 
of the equations 

(ai + a 2 + ni +j)U mJ = (d + a 2 )(l - 5 y o)II m> _,-i 

+ »i(l - 5mo)n m _ w + (7 + l)II m ,, +1 , (73) 
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which satisfy the initial conditions 

n m0 = S m0 , m = 0, 1, • • ■ . (74) 

It follows from (73) and (74), by induction ony, that 

II m> =0, m>j. (75) 

It may be verified, from (72) to (74), that 

Pu = X p m oU. m -ij, Hi + 1 S i < ki, <y < n 2 . (76) 

m — i 

The quantities II mj may be calculated sequentially from the recurrence 
relations (73), with the help of (75) and the initial condition IIoo = 1. 
In order that the representations in (23) and (76) agree fory = n 2 , 
we must have 

*, *, 

£ Pm0n m -,> 2 = £ brSn^Kr, a^S^i-Kr, a 2 )di(K r ) (77) 

m—i r— 

for m + 1 < i < ki. But, from (73), (74), and (93), it follows thatIIo,„ 2 
= Sn 2 (n\, a\ + a 2 ) > 0. Hence, (77) may be used to solve successively 
for /?*,,<), • • • ,p«,+i,o in terms of b r , r = 0, • • •, k\. It remains to satisfy 
the boundary conditions at i = n\ , <y < n 2 — 1 and at j = n 2 ,n\<i 
< k\ , as well as (34) for < i < ni — 1, and the normalization condition 
(71). From (70), these boundary conditions are 

(ai + a 2 + n\ + j)p ni j = o-\Pn^-\j 

+ (ai + a 2 )(l - 5,o)pn„>-i + nxp nx+l j + (y + l)p„,, y+1 , (78) 

for 0<y'< « 2 - 1 and 

Oi(l - fi,- Al ) + a 2 + n x + /i 2 ]p,,„ 2 

= aipi-i,^ + (ai + a 2 )A> 2 -i + «i(l - fcO/X+i.* + n 2 pi,^+i, (79) 

for Aii < i < ^i . 

From (27), (76), and (78), we find, after reduction, that 

X brOn^Kr^KrSn^l + K r , Cll)Sj(-Kr, Cl 2 ) + OlSn^Kr, a } )Sj{-l - K r , do)] 
r-0 

"2 

+ £ c,[r},s ni (l + i)/, a i )Sj(—q,, a 2 ) + aiS ni (Tj/, ai)s y (-l - tj,, a 2 )] 



.-i 



- nx £ Pmon m - n ,-i,y = o, (80) 
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for <> < n 2 - 1. Also, from (23), (76), and (79), it is found that 

02 X fergn^^r, ai)S Wa (-ICr,a2)g|(Kr) . , > 

- (ai + a 2 ) 2 p B «oII n t-»,« a _i = 0, (81) 

m = i 

for /ii + 1 < i < ki . 

From (23) and (27), the boundary condition corresponding to i = ni 
in (79) reduces to 

02 2j b r S ni (Kr, ai)S^(-JC r , flfe/ft^ («r) 



■0 <J»n 2 (Kr) 

- (ai + a 2 ) J ftrS„,(K r , aOs^-^-Kr, a 2 )dn l (K r ) 



n- 2 



- (ai + a 2 ) £ c*8fc,(i|/, ai)s„ 2 _i(-Tj^, a 2 ) = 0. (82) 



This condition may be shown to be redundant, by summing (80) from 
j = to n 2 — 1, (34) from i = to ni — 1, and (81) from i = n\ + 1 to k\ , 
and adding. 

The constants 6 r , r = 0, • • ■ , k\ and c,r, <? = 1, • • • , /i 2 , and the 
probabilities p m0 , m = ni + 1, • • • , ki, are determined by (34), (77), 
(80), and (81) only to within a multiplicative constant, which is 
determined by the normalization condition (71). But, from (70), if we 
sum on j from to k 2 , and on i from «f to k\ , we obtain 

ni £ Po- = ai £ p,- W , »i + 1 < <f < ki. (83) 

Hence, from (23), (27), and (83), with the help of (11), (15), and (98), 
the normalization condition (71) implies that 

£ b r S ni (l + K r , Cl\)0nMr) 



Sno{-Kr, a 2 ) ; . . + Snj-Ul ~ K r , a 2 ) 
<Pn,(K r ) 



"2 



+ X c,s„,(l + tj,, ai)Sn 2 -i(l - t)t, a 2 ) 



+ -Z b r Sn s (Kr, a 1 )Sn 2 (-K rf a 2 )[^ q ,(K r ) ~ 1] ^^T " L < 84 > 

VI. SOME STEADY-STATE QUANTITIES: CASE II 

We now consider the calculation of various steady-state quantities 
of interest. Since no overflow is permitted from the primary queue, R i2 
= in Fig. 1. The loss probabilities Li and L 2 are given by (41), and 
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the probabilities Qi and Q> that demands from the primary and 
secondary sources are queued on arrival are given by (42). Hence, from 
(23) and (27), it follows that (53) to (56) still hold. We emphasize, 
however, that the constants b r , as well as c ( , differ for the two cases. 
Similarly, R 22 is given by (46) and V 2 is given by (48), so that R 22 = 
A 2 <?2, as before, and (62) holds. Also, 

Xi = 2 2 mind", ni)pij = m - 2 J (n, - i)p ijy (85) 

i=o y-o i-o y=o 

and it follows, from (23) and (27), that (64) holds. 

The probability that a demand arriving from the primary source 
overflows (immediately) is 

*1 n 2 -l ti 2 -1 / A, \ 

Ji2= 2 2 Pij= 2 [Pn i j+ 2 Pu)> (86) 

,=n, y'-0 y=0 \ i=n,+l ) 

Hence, with the help of (83), we obtain 

"2-1 *1 k 2 / \ 

Iv= 2 Pn t j+ 2 2 — Pi-Ij-Pij . (87) 

y-0 i-n,+1 y-nj \/li / 

From (23) and (27), with the help of (11), (15), and (98), it is found 
that 



A, 



I\2 = 2 brSn^Kr, ai)s^-i(l - fC r , C^)^, (>Cr) 
r-0 

"2 A, 

+ 2 <vs„,(7j,, ai)s„ 2 _i(l - 7j,, a 2 ) + 2 6 r s„.(K r , ai)s^(-K r , a 2 ) 

f-\ r=0 

.h[a Qi ( Kr) -i]-a q ^Kr)}^l. < 88 > 

Since R i2 = and Ai = ai/x, the first relationship in (52) becomes 

ai(l -Li - I 12 ) =Xi. (89) 

In view of (53), (64), and (88), this relationship provides a useful 
numerical check. 

The mean departure rate from the primary queue to the primary 
servers is 

Rn = n i n 2 2 Pij- (90) 

i = n,+l y— 

From (42) and (83), it follows that R u = XiQi, as expected. Finally, 
the average number of demands in the primary queue is 

^■=22 (i-ni)p tJ . (91) 

i = n, + l y-o 
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Hence, with the help of (11), (15), (23), (83), (113), and (115), it is found 
that 

Vx = bos n <0, Ol K(0, a 2 )^Y'A J-^ ** <0) 



+ — X brSnAKr, ai)S„ 2 (-»C r , Ct 2 ) - . . 



.|n g ,(K r )-(gi + l)+^ 
where A g (£) is as defined in (61). 



7 (, c > | -^Qfc-idBr)!!, (92) 



Oi 



APPENDIX 

We define s,(A, a) by the recurrence relation 

(a + i + X)si(X, a) = a(l - 5,o)s»-i(A, o) + (i + l)s, + i(A, a); 

s (A, a) = 1, (93) 

for i = 0, 1, • • • . Thus s n (X, a) is a polynomial of degree n in both A 
and a, and it may be related to a Poisson-Charlier polynomial. 4,5 
However, we will give here the properties of s„(X, a) which we will 
need. An explicit formula is 1 

* (X ' o) " r 5,r!(i-r)!' 



(94) 



where 

(A) =l, (A) r = A(A + l) ••• (A + r-1), r = 1, 2, ■ • • . (95) 
It was also shown 1 that 

(i + l)s I+ i(A, a) = asi(X, a) + As,(A 4- 1, a) (96) 

and 

s,(A, a) = Si(X + 1, a) - (1 - S,o)s,-i(A + 1, a). (97) 

From (97), it follows that 

£ Si (X, a)=s n (X + l, a), (98) 

«=o 

and, from (96) and (98), we deduce that 

Y (n- i)s t (X, a) = (1 - S n o)s„-i(A + 2, a). (99) 

We now turn our attention to the Chebyshev polynomials of the 
second kind, 3 U, (x). They may be defined by the recurrence relation 

2xUAx) = U, +l (x) + U,-Ax); £/_,(*) = 0, U (x) a 1, (100) 
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for ^« 0, 1, ... . From (10) and (100), it follows that 

(a 2 + n 2 - k)*Ak) = a 2 ¥,+i(ic) + n 2 */-i(/c); 

*_,(*) =0, ¥ (k)«1. (101) 

From (11) and (101), we deduce that 

[a 2 (l - 8jk 2 ) + n 2 - k]4>j(k) = a 2 «fe,-i(ic) + n 2 (l - 8/kJto+iM, (102) 

forj<k 2 . Since 3 



Vt 



IH 



r=0 



it follows that 



-i = |sVVs+*) = s fa 

.°2/ \2Va 2 n 2 / ^=0 \a2 



Hence, from (11), we have 



*"»-£ 



*a-y 



Now 3 



EM*) = 



[<* + ^?^T)' +1 - u - v9^ny +1 ] 



2V*^1 



(103) 



(104) 



(105) 



X f* 1. (106) 



Hence, from (10), we obtain 



*,,(*,) H^ 



*/2 



«, \ [(cr + v^TT)' +1 - (e r - J£=l)' +l ] 



2v^T 



, (107) 



for r = 1, • • • , ki, where e r is given by (20). We note that e r > 1 for 
r = 1, • • • , ki, since K r < for r = 1, • • • , k\. We take the positive 
square root of e 2 r — 1 in (107). From (11), it follows that 



fr(Kr) = 



(e r + Ve?- 1 - y/a 2 /n 2 ) 
2sfIT=l 
(e r - Ver - 1 - -Ja 2 /n 2 ) 

2vS^I 



^( er +V^l) 

a 2 



-i*2 -y 



»2 



1*2 "> 



a 2 



, (108) 



for r = 1, • • ■ , fa. Since K r < for r = 1, • • • , k x , we deduce from (20) 
that 



>■■■ ^>- "—i^WShVS < 109) 




for r = 1, • ■ • , k\. 
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Next, from (14) and (100), it follows that 
(a x + Tii + k)£2/(k) = aifi/+i(ic) + hi£2/-i(k); 

8-iOc) - 0, fioU) ■ 1. (110) 
From (15) and (110), we deduce that 

[aid - «.*,) + ni + kWk) = Oift-i(ic) + »i(l - 8 ikl )8 i+ i( K ), (111) 
for i < k\. Also, from (14) and (103), we obtain 

Ri\ _ T /ai + /ii\ _, //ii 



fl ' (0)=( W %^J=IU'- (112) 



Hence, from (15), we have 

k,-i 



ai 



em = ( - 1 . (ii3) 



Since k\ = rt\ + q\, it follows that 



| (i - Wl )ft(0) = (-YX, (-) , (114) 

!=«! \a\J \n\J 

where A 9 (|) is as defined in (61). Finally, if we multiply (111) by 
(i — ni) and sum on i, we obtain, with the help of (15), 

k £ (i - ni)Bi(K) = ai[J2 9l (»c) - 1] - rtiB, r ,(/c). (115) 
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